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ABSTRACT 


The  project  has  been  successful  In  predicting  decreased  erosion  of  surfaces 
exposed  to  high  heat  fluxes  under  railgun  or  EM  launcher  conditions.  The  signi¬ 
ficant  reduction  has  come  as  a  result  of  the  proper  modeling  of  the  vapor 
shielding  (VS)  effect.  Favorable  comparison  with  experimental  results  supports 
our  conclusion.  A  further  decrease  In  surface  erosion  should  come  as  a  result 
of  an  applied  magnetic  field  which  should  reduce  turbulence  In  the  VS  boundary 
layer.  Preliminary  experimental  measurements  support  this  claim,  but  further 
theoretical  work  must  be  done  which  would  include  2-0,  and  turbulent  effects. 
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MAGNETIC  VAPOR  SHIELD  MECHANISM  FOR  PROTECTION 


OF  MATERIAL  SURFACES  EXPOSED  TO  HIGH  HEAT  LOAD 
Contract  Period  (January  1986  -  September  1989) 

John  G.  Gilligan 

Department  of  Nuclear  Engineering 
North  Carolina  State  University 

INTRODUCTION 

This  final  report  summarizes  theory  and  numerical  modelling  work  accom¬ 
plished  on  our  project  entitled  "Magnetic  Vapor  Shield  Mechanism  for  Protection 
of  Material  Surfaces  Exposed  to  High  Heat  Load"  (DAAL03-86-K-0029).  The  time 
period  of  the  contract  ran  from  January  1986  through  September  of  1989  with 
total  support  of  approximately  $151,000.  A  two-year  contract  renewal  proposal 
(~  $120,000)  has  been  submitted  (April  1989)  and  is  currently  being  evaluated  by 
the  AR0.  Support  for  the  project  has  come  from  both  the  engineering  (David 
Mann)  and  materials  (Andrew  Crowson)  divisions  at  AR0.  A  larger  companion  pro¬ 
ject  for  experimental  verification  and  testing  of  the  magnetic  vapor  shield 
(MVS)  mechanism  has  been  funded  by  the  SDI/IST  through  the  ARO  at  a  level  of 
$300,000  from  July  1987  through  December  1989.  It  is  likely  that  fundinr  of 
that  project  (through  the  SDC)  will  continue  at  a  level  of  $612,000  from  January 
1990  through  December  1993. 

Overview 

Plasma-material  surface  interaction  is  of  vital  importance  to  research  areas 
supported  by  the  Department  of  Defense.  Plasma  etching  and  vapor  deposition 
techniques  are  used  in  the  fabrication  of  microelectronics  devices.  Combustion 
plasmas  are  present  in  propellant  systems  and  engines.  Vaporization  plasmas  are 
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created  when  material  surfaces  are  exposed  to  high  heat  flux  in  gun  barrels, 
plasma  switches,  EM  launchers  (rail guns)  and  electrothermal  guns.  The  goal  of 
our  research  is  to  computationally  model  and  experimentally  verify  the  Magnetic 
Vapor  Shield  (MVS)  concept1”2  as  a  means  of  protecting  material  surfaces  from 
ablation  and  melting  caused  by  extremely  high  heat  loads.  The  basic  vapor 
shield  mechanism  has  been  observed  by  our  group  in  an  electrothermal  gun2 
(SIRENS)  and  by  Keefer  in  an  EM  launched.  The  important  role  of  ablation  in 
rail  gun  operation  has  been  reviewed  by  Parker*. 

The  basic  senario  can  be  visualized  as  exposing  selected  sample  materials 
(insulators,  conductors,  semiconductors)  to  a  plasma  source  equivalent  to  the 
high  heat  fluxes  experienced  in  many  applications  such  as  gun  barrels,  switches 
and  railguns.  A  simple  vapor  shield  (See  Fig.  1)  resulting  from  ablated  surface 
material  is  fairly  common  when  high  heat  or  particle  fluxes  are  present. 

However,  the  idea  of  using  magnetic  fields  to  actively  control  transport  coeffi¬ 
cients  and  suppress  turbulence  under  these  extreme  conditions  is  relatively 
unexplored. 1-2  The  critical  questions  to  be  addressed  are  whether  reducing  the 
vapor  thermal  conductivity  (dominated  by  electron  collisions)  will  decrease  the 
total  heat  load  on  the  material  surface  and  whether  magnetic  fields  can  be  made 
sufficiently  strong  so  as  to  reduce  cross-field  turbulence.  Calculational  evi¬ 
dence  obtained  thus  far  suggests  that  both  questions  can  be  answered  in  the 
affirmative. 

Of  special  interest  to  the  ARO  has  been  the  ablation  of  internal  insulators 
and  electrodes  exposed  to  the  radiation  from  the  accelerating  plasmas  in  railguns. 
These  devices  have  strong  self-induced  magnetic  fields.  The  MVS  concept  pre¬ 
dicts  that  if  sufficient  magnetic  field  is  present  at  the  material  surface, 
then  the  ablated  plasma  transport  and  covection  coefficients  can  be  drastically 
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reduced  resulting  in  less  energy  being  transferred  to  the  surface.  In  turn,  the 
energy  is  re-directed  into  axial  convective  particle  transport  or  photon 
radiation  transport  away  from  the  localized  area.  It  is  obvious  that  the  viabi¬ 
lity  of  the  MVS  concept  is  dependent  on  the  delicate  balance  between  the  energy 
flows,  electromagnetic  radiation,  particle  convection,  and  conduction.  Hence, 
for  a  computational  analysis,  the  coupled  energy,  momentum  and  particle  balance 
equations  must  be  solved  simultaneously  and  self-consistently. 

Goals  of  Project 


•  To  understand  and  control  the  basic  energy  transport  phenomena  in  the  VS 
at  a  material  surface  under  high  heat  flux. 

•  To  explore  the  effect  of  a  strong  magnetic  field  in  decreasing  surface 
erosion  under  high  heat  flux. 

•  To  model  the  MVS  mechanism  and  to  compare  with  experimental  results. 

Principal  Project  Accomplishments  (1986-1989) 

•  Described  basic  physical  mechanisms  of  the  vapor  shield  (VS)  process  for 
protection  of  surface  conditions  including  necessary  conditions. 

•  Described  basic  physical  mechanisms  to  enhance  the  VS  effect  with  strong 
magnetic  field  -  Magnetic  Vapor  Shield  (MVS). 

•  Developed  1-D,  time  dependent  heat  transfer  and  phase  change  code 
(MELTIT).  ‘ 

•  Developed  1-D,  time  dependent,  MHD,  radiation  transport  code  (MAGFIRE)  to 
analyze  VS  and  MVS  effects. 

•  Developed  0-D,  time  dependent,  energy  and  particle  balance  code  (ZEUS)  to 
compare  theoretical  results  with  experimental  data. 

•  Demonstrated  that  radiation  transport  is  one  of  the  most  important  energy 
transfer  mechanisms  in  the  VS.  Further,  that  flux-limited  diffusion 
theory  is  adequate  to  describe  the  transport. 

•  Demonstrated  with  numerical  simulation  that  the  VS  is  much  more  effective 
in  reducing  erosion  than  previously  throught  (-  10%  energy  transmission 
through  the  VS  as  opposed  to  -  50%). 

•  Demonstrated  that  certain  materials  are  resistant  to  high  heat  flux 
damage,  namely  some  ceramics  and  refractory  metals.  Some  materials  have 
better  VS  properties  than  others. 
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•  Simulation  of  the  MVS  mechanism  was  inconclusive  in  1-D  although  experi¬ 
ment  has  shown  a  decrease  in  erosion  with  increasing  magnetic  field. 
Simulation  with  2-0  is  required. 

•  Comparison  of  numerical  simulation  and  experiment  has  been  excellent  for 
prediction  of  erosion  depths. 

•  Eleven  conference  presentations,  nine  refereed  publications  and  one 
invited  talk  have  been  given  on  the  vapor  shield  project  (See  listing  of 
publications  at  end  of  report).  Dr.  J.  Gilligan  was  awarded  the  NCSU 
College  of  Engineering  Alcoa  Research  Prize  for  work  on  the  project  over 
the  past  three  years. 

•  Two  masters  theses  and  one  Ph.D.  thesis  have  been  produced.  Two  current 
graduate  students  are  partially  or  fully  supported  by  the  project. 

M.S.,  1986,  M.  Holliday,  "Analysis  of  Surface  Material  Ablation  and 

Melting  Using  the  Heat  Balance  Integral  Technique"  (J.  Gilligan  and  M. 
Doster,  Advisors). 

M.S.,  1990,  R.  Rustad,  "Melting  and  Ablation  of  the  First  Wall  During  a 
Plasma  Disruption"  (J.  Gilligan,  Advisor). 

Ph.D.,  1990,  D.  Hahn,  "Energy  Transport  Through  a  Plasma  Boundary  Layer" 
(J.  Gilligan,  Advisor). 


PROJECT  RESULTS  AND  DISCUSSION 

Before  the  beginning  of  this  project  the  concept  of  vapor  shielding  was 
established  in  other  fields^-7  put  few  detailed  calculations  had  been  performed 
to  analyze  its  effectiveness.  Early  estimates  for  the  total  energy  fraction 
transmitted  through  the  vapor  shield,  f,  were  in  the  range  of  50%. 6  However, 
these  results  did  not  necessarily  agree  with  experimental  measurements  for  laser 
light  incident  on  surfaces  which  were  much  lower  (10-20%)7.  The  first  task  of 
this  project  was  to  verify  that  under  the  typical  conditions  experienced  in 
railguns  (heat  flux,  duration,  photon  spectrum,  etc.)  that  the  VS  would  be  opti¬ 
cally  thick  enough  to  absorb  the  Incident  photon  radiation.  With  simple 
estimates  for  the  vapor  production  rate  and  detailed  calculations  of  photon 
radiation  opacities,  we  were  able  to  determine  that  there  was  more  than  suf¬ 
ficient  attenuation  of  photons  in  the  VS  to  reduce  the  incident  heat  flux.1  In 
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addition,  it  was  found  that  magnetic  fields,  if  sufficiently  strong,  could 
reduce  cross  field  thermal  transport  dramatically  if  turbulence  could  be 
suppressed  as  had  already  t>een  noted  for  non-ablative  boundary  layers.8 

In  order  to  confirm  the  initial  results,  several  codes  were  developed  to 
properly  analyze  the  high  heat  flux,  plasma  boundary  layer  at  a  surface.  The 
codes  and  their  interrelationships  are  briefly  illustrated  in  Figure  2.  Codes 
that  have  been  completed  and  tested  are  indicated  in  solid  blocks.  Codes  under 
development  are  described  in  blocks  with  dashed  lines.  Interfaces  with  experi¬ 
mental  results  are  also  indicated.  Individual  code  results  will  now  be 
discussed  for  the  vapor  shield  effect. 

MELTIT:  1-D,  Time  Dependent,  Phase  Change  and  Heat  Transfer 

In  order  to  calculate  the  amount  of  surface  material  vaporized  and/or  melted 
under  high  heat  flux,  the  1-D,  time  dependent  code  MELTIT  was  written. 9  The 
heat  balance  integral  technique  was  used  to  model  the  temperature  distribution 
throughout  the  surface  material.  Phase  change  and  resulting  moving  boundaries 
were  explicitly  calculated.  In  addition,  the  effects  of  molten  material  ejec¬ 
tion  and  vapor  redeposition  were  included.  However  vapor  shielding  was  not 
directly  modeled  in  MELTIT.  Figure  3  shows  typical  results  generated  with  the 
code  giving  the  melt  front  and  left  boundary  position  as  a  function  of  time  for 
an  incident  energy  density  of  5x106  J/m2  on  a  0.5  cm  slab  of  SS  316.  The 
effects  of  melted  material  ejection  can  be  seen  where  25%  of  the  molten  material 
formed  is  assumed  to  be  removed.  As  expected,  ejection  of  molten  material  not 
only  increases  the  erosion  rate  at  the  left  boundary,  but  increased  the  melt 
front  velocity  since  there  is  less  resistance  to  heat  transfer.  The  ablation 
rate  calculated  from  MELTIT  is  an  important  variable  that  determines  the  VS 
transmission.  Different  materials  were  found  to  experience  more  ablation  for  a 
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given  incident  energy  fluence,  if  they  had  low  thermal  conductivity  and  low 
sublimation  energies  (typical  of  insulators  and  some  ceramics). 

HAGFIRE:  1-D,  Time  Dependent,  MHO,  Radiation  Transport 

To  analyze  the  VS  mechanism  requires  simultaneous  solution  of  models  of: 

•  Heat  conduction  and  phase  change  in  the  material  surface, 

•  MHO  transport  of  the  shielding  vapor,  including  electric  and  magnetic 
field  effects, 

•  radiation  transport  in  the  shielding  vapor, 

•  incident  radiation  energy  and  momentum  deposition, 

•  equation-of-state  data  for  the  shielding  vapor. 

The  code  MELT1T,  previously  described,  satisfied  the  first  modeling  require¬ 
ment.  For  the  transport  of  the  ablated  vapor  we  employed  a  1-D,  time  dependent 
Lagrangian  MHD  code,  ZPINCHlO,  which  incorporates  multi-group  radiation  trans¬ 
port.  ZPINCH  was  adapted  to  account  for  vapor  creation  at  the  boundary  surface, 
low  temperature  collisions,  a  parallel  external  magnetic  field  and  high  magnetic 
field  effects  on  transport  coefficients,  and  was  renamed  MAGFIRE.il 

For  the  following  results,  the  heat  flux  was  modeled  as  a  blackbody  spectrum 
(divided  into  twenty  energy  groups)  with  a  constant  intensity  given  by 

Os  ■  oTs. 

where  o  is  the  Stefan-Boltzmann  constant  and  Ts  is  the  source  temperature. 

Values  of  T$  in  the  range  2-5  eV  were  chosen  which  give  heat  fluxes  of 
(1.6-64.0)  X  1010  w/m2  that  are  typical  of  railguns.  The  heat  flux  is  incident 
on  a  material  surface  with  a  characteristic  heat  of  vaporization. 

Shown  in  Figs.  4  and  5  are  density  and  temperature  profiles  for  the  vapor 
shield  as  a  function  of  distance  from  the  ablating  iron  surface  (semi-infinite 
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slab  geometry)  at  times  of  1,  10,  and  100  us.  The  incident  heat  flux  is 
8.3  X  1010  W/m2  corresponding  to  a  black  body  with  Ts  ■  3  eV.  The  density 
always  remains  high  immediately  adjacent  to  the  ablating  surface  and  then  falls 
off  as  the  outer  vapor  shield  layers  expand.  The  temperature  of  the  shield 
vapor  is  low  near  the  surface  but  then  increases  away  from  the  surface  due  to 
absorption  of  the  incoming  radiation.  It  should  also  be  noted  that  the  average 
ionization  state  was  found  to  change  -  10-2  at  the  surface  to  over  2.0  in  the 
high  temperature  zones. 

The  ablation  erosion  thickness  for  various  incident  radiant  fluxes  and 
materials  is  shown  in  Table  I  for  10  us  pulses.  Also  given  is  the  percent  of 
transmitted  energy  f  through  the  vapor  shield  to  the  material  surface.  Previous 
estimates  of  f  not  using  radiation  transport  techniques  or  MHD  analysis  of  the 
transmitted  energy  were  in  the  50%  range. 6  Our  results,  where  the  transmitted 
energy  ranges  from  21%  down  to  3%,  presented  a  much  different  picture  of  vapor 
shield  effectiveness.  Thus,  the  predicted  erosion  thickness  can  be  up  to  an 
order  of  magnitude  less  than  that  estimated  with  a  crude  vapor  shield  model. 
Typical  ablation  thicknesses  are  on  the  order  of  microns  for  various  materials 
in  Table  I.  For  higher  heat  fluxes  the  transmitted  energy  fraction  decreases 
because  more  vapor  is  produced  and  it  provides  a  "thicker"  shield.  No  clear-cut 
erosion  resistant  material  is  evident  in  Table  I.  Carbon  has  a  relatively  high 
heat  of  vaporization  (7.4  eV)  compared  to  iron  (4.3  eV)  and  copper  (3.5  eV). 
fhus,  the  carbon  vapor  shield  is  somewhat  less  effective  for  a  given  heat  flux 
since  less  vapor  is  produced. 

The  effect  of  ablation  rate  on  transmitted  energy  through  the  vapor  shield 
is  shown  in  Fig.  6.  The  source  spectrum  is  fixed  for  a  3  eV  blackbody  but  the 
intensity  is  increased  or  decreased  to  give  different  incident  radiation  fluen- 
ces  for  a  fixed  time  of  10  ps.  The  actual  fluence  for  a  3  eV  blackbody  over  10 
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ys  Is  8.3  105  j/m^.  As  can  be  seen,  the  larger  fluences  give  lower  transmission 
fractions  through  the  vapor  shield.  This  is  a  result  of  the  larger  amount  of 
vapor  produced  which  makes  the  vapor  shield  "thicker".  It  appears  there  may  be 
saturation  of  the  transmitted  energy  at  higher  fluences.  Likewise,  there  may  be 
a  cutoff  below  which  no  vapor  shielding  takes  place.  However,  the  validity  of 
the  radiation  transport  model  that  uses  the  diffusion  approximation  is  in 
question  for  the  low  fluence  cases  below  104  J/m2.  In  addition,  less  vapor  is 
actually  produced  at  the  low  fluences  than  in  our  vapor  source  model  since 
melting  would  become  more  important. 

Results  from  our  detailed  simulations  have  revealed  that  photon  radiation 
transport  (absorption  and  re-emission)  is  the  dominant  means  of  energy  transport 
through  the  vapor  shield.  Shown  in  Figure  7  is  the  energy  fluence  reaching 
the  material  surface  as  a  function  of  time.  Radiative  energy  is  clearly  domi¬ 
nant.  However,  these  simulations  only  include  simple  heat  conduction  through 
the  shield  which  gives  rise  to  a  very  small  contribution  to  the  total  energy 
fluence.  Observations  by  other  researchers  indicate  that  convective  turbulence 
may  be  as  important  as  radiation  in  the  transfer  of  energy  to  the  surface. 12 
Some  of  our  preliminary  simulations  with  "enhanced"  convective  transport  have 
indeed  indicated  that  this  may  be  the  case. 

Since  radiation  transport  plays  such  an  important  role  in  the  vapor  shield, 
we  have  re-examined  some  of  the  basic  principles  used  in  the  radiation  transport 
modeling.  In  MAGFIRE  the  radiation  diffusion  approximation  is  used.  The  basic 
assumption  of  the  diffusion  approximation  is  that  the  radiation  field  is  isotro¬ 
pic.  However  this  is  valid  for  an  optically  thick  medium  only.  Near  the 
surface  of  an  optically  thick  medium,  or  for  optically  thin  systems,  the 

radiation  diffusion  equation  overestimates  the  radiation  transport.  Thus,  the 
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diffusion  approximation  overestimates  the  radiation  flux  at  the  wall  surface  due 
to  the  boundary  effect,  regardless  of  the  optical  property  of  the  vapor  shield 
as  a  whole.  One  way  to  overcome  the  problem  of  overestimation  is  to  use  a  flux 
limiter.  The  idea  is  to  alter  the  diffusion-theory  formula  for  the  flux  in  such 
a  way  as  to  yield  the  standard  result  in  the  high-opacity  limit,  while  simu¬ 
lating  free  streaming  in  transparent  regions. 

A  formal  solution  for  the  radiative  transfer  equation  was  derived  in  order 
to  compare  with  the  present  numerical  solution  which  uses  the  diffusion  approxi¬ 
mation.  One  can  obtain  an  analytic  solution  of  the  radiative  transfer  equation 
for  the  steady  state  case  where  the  temperature  and  density  distributions  and 
also  the  radiation  field  are  independent  of  time.  It  is  assumed  that  the 
radiation  field  is  weakly  dependent  on  angle  to  derive  the  multigroup  form  of 
the  radiative  transfer  equation.  By  integrating  the  solution  Ig(n)  over  all 
solid  angles,  one  can  calculate  the  radiant  energy  flux  Sg  at  the  wall  surface. 
In  order  to  compare  the  equivalent  quantities,  our  one  dimensional,  time 
dependent,  MHD  code  MAGFIRE  is  run  first  until  quasi -steady  state  is  reached, 
and  then  by  using. the  profiles  of  density  and  temperature  calaulated,  the 
analytic  solution  for  the  steady  state  multigroup  radiative  transfer  equation  is 
obtained.  The  analytic  solution  consists  of  two  parts,  one  for  the  attenuation 
of  the  external  source  radiation,  and  the  other  for  the  radiation  from  the 
heated  vapor  shield  plasma. 

As  shown  in  Table  II,  the  diffusion  model  with  an  isotropic  flux  limiter 
(c/4,  c=  speed  of  light)  underestimates  the  radiation  flux  compared  to  the 
analytic  solution.  The  maximum  anisotropy  value  of  c  overestimates  fluxes  and 
further  overestimation  would  be  obtained  without  a  flux  limiter  at  all.  By 
using  the  flux  limiter  value  of  (c/2),  there  is  a  relatively  good  agreement 
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between  MAGFIRE  and  the  analytic  solutions.  This  is  due  to  the  fact  that  the 
radiation  field  is  somewhat  anisotropic  and  has  a  streaming  behavior  inside  the 
vapor  shield  plasma,  especially  near  the  wall.  Figs.  8  and  9  show  the  frequency 
group  dependent  radiation  fluxes  at  10  psec  for  the  flux  limiters  of  (c/4)  and 
(c/2)  respectively. 

Since  MAGFIRE  gives  the  radiation  flux  as  a  result  of  diffusive  processes, 
it  is  impossible  to  distinguish  the  attenuated  external  radiation  source  from 
the  radiation  emitted  and  attenuated  by  the  vapor  shield  plasma.  However,  ana¬ 
lytic  solutions  show  that  the  attenuated  source  flux  decreases  with  time  due  to 
the  increased  absorption  by  the  vapor  shield  plasma,  and  radiation  from  the 
vapor  shield  increases  with  time  due  to  the  increased  temperature  of  the  vapor. 

The  flux  limiter  of  (c/2)  was  chosen  for  the  diffusion  model,  but  the  net 
effect,  that  is,  the  amount  of  erosion,  is  insensitive  to  the  choice  of  the  flux 
limiter  values.  Transmission  factors  for  the  case  of  the  iron  wall  with  3  eV 
blackbody  radiation  of  10  psec  pulse  are  8.9,  10.8,  and  12.3%,  for  the  flux 
limiter  values  of  (c/4),  and  (c/2),  and  (c),  respecti vely.  The  small  differen¬ 
ces  in  transmission  factors  are  due  to  the  self  limiting  mechanism  of  the  vapor 
shield,  that  is,  the  larger  energy  flux  causes  more  ablation,  more  vapor  shield, 
and  then  a  smaller  flux. 

ZEUS:  0-D,  Time  Dependent  Electrothermal  Gun  Simulation 

Since  the  results  from  the  MAGFIRE  code  were  not  directly  applicable  to 
experimental  observation,  the  ZEUS  code  was  written  to  incorporate  MAGFIRE 
results  (i.e.  predictions  for  f)  in  the  simulation  of  the  breech  (source)  region 
of  our  electrothermal  launcher  SIRENS.  ZEUS^  is  capable  of  predicting  erosion 
rates  of  the  breech  region  surfaces  as  well  as  plasma  temperatures,  densities. 
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pressures,  ionization  states,  etc.  and  thus  becomes  the  direct  link  to  experi¬ 
mental  observations.  ZEUS  describes  the  physical  behavior  of  ablation 
controlled  devices  using  the  vapor  shield  concept.  The  ZEUS  code  has 
demonstrated  good  agreement  with  previous  steady-state  ablation  plasma  models 
and  experimental  results  and  also  yields  very  reliable  predictions  for  ablation 
depth  in  SIRENS. 

The  0-D,  time  dependent  code  ZEUS  includes  models  for  the  energy  transport, 
particle  transport,  resistivity  and  equations-of-state  under  local  thermodynamic 
equilibrium.  The  degree  of  ionization  for  the  plasma  is  computed  from  a 
modified  Saha  equation  based  on  arc  temperature  and  density.  The  input  parame¬ 
ters  to  the  code  include  the  initial  density,  temperature,  charge  state  of  the 
plasma  and  the  time  history  of  the  discharge  current.  The  code  has  the  addi¬ 
tional  feature  of  inputting  the  energy  transmission  factor  f  through  the  vapor 
shield  region  at  the  material  surface  as  predicted  by  MAGFIRE.  Special  models 
also  have  been  developed  for  non-ideal  plasma  effects  and  the  calculation  of 
ablated  material  enthalpy. 

A  set  of  measurements  at  different  discharge  energies  was  carried  out  to 
study  ablation  of  the  Lexan  insulator  breech  in  SIRENS  (Fig.  10). 

The  geometric  parameters  for  the  breech  region  of  SIRENS  are: 

Diameter  *  4mm 

Length  *  12  cm 

A  typical  pulse  discharge  current  profile  as  a  function  of  time;  Imax  =  15  kA  at 
t  *  8  us,  pulse  length  100  ys;  was  used  as  the  input  to  the  code.  The  current 
profile  of  this  reference  shot  was  scaled  to  obtain  the  maximum  current  of  dif¬ 
ferent  shots  so  as  to  compare  the  theoretical  and  experimental  results  at 
different  plasma  variables  which  include  density,  temperature,  enthalpy  and 
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charge  state.  The  heat  flux  to  the  insulator  surface  is  determined  from  the 
plasma  temperature  by  assuming  that  a  fraction  of  the  blackbody  radiation,  f, 
impinges  on  the  wall.  The  variable,  f,  was  varied  from  0.2  to  0.01  and  the 
total  ablation  depth  from  the  code  was  compared  to  that  from  the  experiment. 
Corroboration  between  the  experiment  and  theoretical  simulation  of  the 
electrothermal  gun  (SIRENS)  was  achieved  for  f  ->  0.2  at  lower  heat  fluxes  and 
f  -  0.05  at  higher  heat  fluxes;  as  shown  in  Fig.  11.  This  corresponds  to  the 
trends  for  f  predicted  by  MAGFIRE. 

The  critical  parameter  involved  in  the  ablation  process  is  the  fraction  of 
blackbody  radiation  transmitted  through  the  vapor  shield,  f.  Fig.  11  shows  that 
the  ablation  depth  increases  with  an  increase  in  the  fraction  of  blackbody 
radiation  through  the  vapor  shield.  In  the  case  of  the  fraction  of  radiation 
transmitted  through  the  vapor  shield  a  higher  flux  of  photons  reaching  the 
surface  of  an  insulator,  would  produce  a  higher  amount  of  ablation.  Since 
the  heat  flux  to  the  wall  is  proportional  to  the  transmission  factor,  f,  any 
increase  in  f  implies  a  corresponding  increase  in  ablation.  The  input  joule 
heating  energy  is  primarily  distributed  between  radiation  energy  and  internal 
energy  of  the  plasma.  Fig.  J1  indicates  that  the  variation  in  joule  heating 
below  about  2000  J  yields  a  proportional  change  in  radiation  energy.  However, 
above  2000  J  the  internal  energy  of  the  plasma  plays  a  dominant  role  in  the 
energy  balance  and  the  heat  flux  does  not  scale  directly  with  input  energy.  A 
fixed  value  of  the  heat  of  sublimation,  Hjub.  for  Lexan  was  adopted  since  the 
calculation  of  that  parameter  is  based  on  similar  principles  as  that  for  PTFE 
(Teflon). 

Fig.  12  shows  typical  predictions  of  the  average  plasma  density  and  tem¬ 
perature  as  a  function  of  time,  for  two  different  values  of  f(f«0.1  and  f*0.01). 
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Predictions  from  the  code,  ZEUS,  indicates  the  following  trend  with  respect  to 
the  fraction  of  radiation  transmitted  through  the  vapor  shield  to  the  wall  sur¬ 
face.  A  higher  value  of  f  produces  a  larger  amount  of  energy  loss  in  the  form 
of  radiation.  The  incident  radiation  on  the  wall  then  produces  ablation  of  the 
wall  material.  The  cold  particles  from  ablation  entering  into  the  plasma  cools 
it  down  further.  The  higher  radiation  energy  loss  combined  with  the  higher 
ablation  rate  gives  rise  to  a  lower  temperature  and  higher  density  with  f*0.1 
relative  to  f*0.01. 

The  temperature  of  the  bulk  plasma  rises  rapidly  as  the  current  increases 
due  to  ohmic  heating  and  reaches  a  peak  at  about  the  same  time  the  current  does. 
At  this  point  of  the  discharge  the  electric  power  is  large  compared  to  the 
radiation  power.  The  discharge  does  not  dissipate  the  supplied  power  and  heats 
up  rapidly.  The  above  mechanism  may  be  considered  to  be  adiabatic.  The  plasma 
density  also  starts  to  increase  rapidly  due  to  a  rise  in  the  temperature. 
However,  due  to  an  increase  in  the  pressure  within  the  channel,  the  particles 
are  transported  out  in  the  axial  direction.  The  plasma  density  reaches  a  peak 
value  at  about  10  ys  later  from  the  peak  plasma  temperature  due  to  the  competing 
effects  of  particle  removel  by  axial  transport  and  ablation.  After  about  50  ys 
the  temperature  reaches  a  quasi-stationary  state,  when  the  temperature  of  the 
discharge  remains  constant  with  time.  This  occurs  due  to  the  radiation  power 
being  of  the  same  order  of  magnitude  as  the  electric  power.  However,  the  rate 
of  particle  loss  is  still  higher  than  the  rate  of  ablation  and  so  the  density 
keeps  decreasing. 

Non-Ideal  Plasma  Effects 

For  the  new  regime  of  high  density,  low  temperature  plasmas  encountered  for 
our  applications,  several  "non-ideal"  effects  were  addressed;  namely,  a  decrease 
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in  the  effective  ionization  potential  of  the  atoms  in  the  vapor  shield  and  the 
effect  on  transport  coefficients  such  as  resistivity,  thermal  conductivity,  etc. 
The  former  phenomena  will  increase  the  ionization  state  of  the  vapor  and  will 
give  rise  to  a  stronger  magnetic  field  induced  current  and  shielding.  The 
latter  mechanism  has  been  analyzed  by  several  researchers1®*1®.  We  have 
included  these  formulations  in  ZEUS  to  find  that  predicted  ablation  depths  are 
reduced  by  about  10%  as  shown  in  Figure  13.  However,  a  more  detailed  analysis 
of  the  effect  on  equation-of-state  and  opacity  data  is  currently  being  investi¬ 
gated. 


CONCLUDING  REMARKS 

The  project  has  been  successful  in  predicting  decreased  erosion  of  surfaces 
exposed  to  high  heat  fluxes  under  railgun  or  EM  launcher  conditions.  The  signi¬ 
ficant  reduction  has  come  as  a  result  of  the  proper  modeling  of  the  vapor 
shielding  (VS)  effect.  Favorable  comparison  with  experimental  results  supports 
our  conclusion.  A  further  decrease  in  surface  erosion  should  come  as  a  result 
of  an  applied  magnetic  field  which  should  reduce  turbulence  in  the  VS  boundary 
layer.  Preliminary  experimental  measurements  support  this  claim,  but  further 
theoretical  work  must  be  done  which  would  include  2-D,  and  turbulent  effects. 
Numerous  researchers  have  reported  reduction  in  heat  transfer  coefficients  under 
somewhat  different  conditions  (i.e.  lower  Reynolds  numbers,  no  plasma  effects, 
etc.).  The  key  question  that  must  be  addressed  in  future  work  is  the  relative 
role  of  radiation  transport  (reported  here)  and  turbulent  convective  transport 
in  the  overall  energy  transport  through  the  VS  boundary  layer. 


V 
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Table  I  Summary  of  Parametric  Studies  for  Vapor  Shielding  Effect 


Source 

Blackbody  Temp.,  eV 

Heat  Flux,  W /m2 

Pulse  Length,  /xsec 

Energy  Fluence,  J  /m2 

2 

1.6xl010 

10 

1.6xl05 

3 

8.3xl010 

10 

8.3x1 05 

5 

6.4xlOn 

10 

6.4x10® 

Iron 

Energy  reaching  wall,  J  /m2 

3.8xl04 

9.0xl04 

2.6x10s 

I 

Transmission  factor  f,  % 

23.0 

10.8 

4.1 

1 

Erosion  Thickness,  pm 

0.7 

1.5 

4.5 

Graphite  : 

Energy  reaching  wall,  J/m2 

7.8xl04 

1.7xl05 

4.5x10s 

1 

Transmission  factor  f,  % 

47.5 

20.5 

7.1 

Erosion  Thickness,  pm 

1.1 

2.4 

6.4 

Copper 
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Radiation  Flux  (W/cmJ)  at  Iron  Wall  Surface  Subjected  to  3  eV 
Blackbody  Radiation  after  10  /*sec  Pulse 


Flux  Limit 
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a)  Attenuated  3  eV  Blackbody  Radiation 

b)  Radiation  from  Vapor  Shield  Plasma 
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Figure  2.  Code  Interelationahi 
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Plasma  Temperature  Distribution  for  the  case  of  an  Iron  Wall  subjected 


Ablation  Velocity  of  the  Iron  Wall  subjected  to  3  eV  Blackbody  Radiation 
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Fig.  7  Energy  in  Vapor  Shield  for  an  Iron  Wall  with  3  eV 
Blackbody  Radiation 
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Fig.  8  Radiation  Flux  at  Iron  Wall  Surface  Subjected  to 
a  3  eV  Blackbody  Rodiation  (  10  pta,  Flux  Limit  —  (c/4)  U  ) 
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Fig.  9  Rodiation  Flux  at  Iron  Wall  Surface  Subjected  to 
a  3  eV  Blackbody  Radiotion  (  10  Flux  Limit  *  (c/2)  Ug  ) 


Fig. 10  Conceptual  Design  of  the  SIRENS  Experiment 
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Fig.  11  Ablation  Depth  of  Lexan  Insulator  in  SIRENS. 
Symbols  show  the  experimental  data  and  solid  lines  are 
the  result  of  ZEUS  calculations. 
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Fig.  13  Ablation  Depth  as  Predicted  by  ZEUS  for  Lexan  with  f~0.05 


